Strong renormalization of the electronic band gap due to lattice polarization in the 

GW formalism 
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The self-consistent GW band gaps are known to be significantly overestimated. We show that this 
overestimation is, to a large extent, due to the neglect of the contribution of the lattice polarization 
to the screening of the electron-electron interaction. To solve this problem, we derive within the GW 
formalism a generalized plasmon-pole model that accounts for lattice polarization. The resulting 
GW self-energy is used to calculate the band structures of a set of binary semiconductors and 
insulators. The lattice contribution always decreases the band gap. The shrinkage increases with 
the size of the longitudinal-transverse optical splitting and it can represent more than 15% of the 
band gap in highly polar compounds, reducing the band-gap percentage error by a factor of three. 



For the past decades we have witnessed a steady in- 
crease of the accuracy of first principles electronic band 
structure calculations. At the lowest level we have den- 
sity functional theory (DFT), normally with a local- 
density (LDA) or a generalized-gradient approximation 
(GGA) to the exchange-correlation functional. Unfor- 
tunately, the Kohn-Sham calculations strongly underes- 
timates quasi-particle band gaps, often by more than a 
factor of two. At the next step, we have the GW approx- 
imation of many-body perturbation theory i^r— For many 
years, the standard practice was to start from a DFT 
calculation, and to evaluate perturbatively the GW en- 
ergy corrections to the Kohn-Sham band structure. This 
procedure, which we will refer to as GM^OLDA, is jus- 
tified only when the departure wave functions and band 
structure are already close to the quasi-particle ones. It 
yields a very good agreement with experimental photoe- 
mission data, especially for standard sp materials. How- 
ever, it is well known that GVF@LDA fails for many crys- 
tals that have d electrons participating in the states close 
to the band gap4i"— To solve this problem one can per- 
form restricted self-consistent (sc) GW , using for exam- 
ple the quasi-particle (QP) scGVF— i^ method or pertur- 
bative GW on top of sc Coulomb hole plus screened ex- 
change4riii£ The scGW^ techniques have the advantage 
of being independent of the starting point and of giving 
accurate results for systems that are not treatable using 
standard approximations, at the price of a larger com- 
putational cost. We should emphasize that these tech- 
niques allow for a precise treatment of materials that were 
considered out-of-reach for ah initio calculations only 10 
years ago. 

Unfortunately, self-consistent GW calculations give 
too large band gaps, especially for polar materials, the 
simplest example being probably LiFi^ The overestima- 
tion of the band gap within QPscGVF was also discussed 
in literature for IITV compounds^^ Zn-IV-N2r^ and 
transition metal oxides j^i^iii^ The causes for these er- 
rors were attributed to the lack of vertex corrections^^ 
or to the underscreening by the random-phase approx- 
imatior^ii (RPA), and in both cases a simple empirical 



correction was adopted to improve the theoretical calcu- 
lations (such as scaling the final self-energy Sxc by about 
80%, obtaining the so-called 0.8E approximatior^i^) . 

The study of electronic correlation and its effects on 
the band-structure of solids has a very long tradition. 
A much less studied subject is, instead, the influence of 
phonons on the quasi-particle spectrum. Only now, in 
fact, the high level of accuracy achieved by electronic 
structure calculations makes the phonon contributions 
larger than the theoretical precision. 

There exist two contributions to the band gap com- 
ing from phonons: i) The electron-phonon coupling, that 
can be handled according to the Allen-Heine-Cardona 
theoryi^ using second-order perturbation theory within 
the harmonic and adiabatic approximations. One obtains 
two terms of the same order to be included in the GW 
self-energy, known as the Fan and Debye- Waller terms, 
ii) The phonon contribution to the frequency-dependent 
dielectric function, that can be a sizable component of 
the total electron screening in polar compounds. 

The electron-phonon coupling terms are responsible for 
the strong reduction of the band gap of diamond. In fact, 
in that case the zero-point contribution to the band gap 
was proved by Giustino et ali^ to be 0.6 eV. They ob- 
served that the example of diamond is an extreme case, 
while for most semiconductors and insulators the zero- 
point renormalization is as small as 10-50 me V— and can 
be safely neglected. The Fan and Debey- Waller terms 
have also been used to determine the temperature de- 
pendence of the band gap of semiconductors J^ 

In this work we want to focus on the second, usually 
overlooked, contribution of phononic nature: the screen- 
ing due to the coupling of the electro-magnetic field and 
the oscillating dipoles i^'^'d^ Of course, this contribution 
exists in polar materials only (and thus it is nonexis- 
tent in diamond). We should emphasize that by "polar" 
we mean any material that has non-negligible Born effec- 
tive charges, leading to a measurable longitudinal-optical 
(LO) and transverse-optical (TO) phonon splitting. Most 
semiconducting and insulating materials known actually 
fit into this category. Indeed, polar materials are char- 



actcrizcd by LO and TO infrarcd-activc phonons, whose 
excitation induces macroscopic electric fields in the crys- 
tal.— These fields contribute, together with the electronic 
screening, to the total screening of the Coulomb poten- 
tial, and their strength is proportional to the bond ion- 
icity. 

Already in 2005 Bechstedt et alt^ observed that the 
lattice polarization affects significantly the dressing of the 
quasi-particles in ionic materials, inducing large correc- 
tions to the band gap, that they calculated using a model 
electronic dielectric function and a static approximation 
for the lattice contribution to the screening. Unfortu- 
nately, it was already observed in Ref. |l3 that that model 
always overestimates the band-gap correction. In this 
Letter, we take a step further and develop a fully ab ini- 
tio framework to include the lattice contributions to the 
screening within the GW scheme. The basic ingredients 
are, besides the standard RPA screening of GW , the LO 
and TO frequencies of the infrared active phonon modes. 
We then calculate the electronic band structure of a set of 
polar materials applying this new scheme, and compare 
the results with experimental data and with many-body 
results including only electronic correlations. 

The lattice polarization plays an important role in the 
determination of the total screening of polar compounds, 
as it is evident from the difference between Eoo, the elec- 
tronic dielectric function at zero frequency, and Cs , the to- 
tal static dielectric function including contributions from 
both electrons and infrared active phonons. The ratio 
between these quantities is related to the phonon fre- 
quencies at the center of the Brillouin zone, wlo and 
wto, through the Lyddane-Sachs- Teller relation: 
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(1) 



The standard treatment of optical phonons in crystals 
using the dynamical matrix implicitly assumes that the 
interatomic interactions are instantaneous. In the case of 
polar crystals, the long-range nature of the Coulomb in- 
teraction requires a proper account of retardation effects. 
It is easy to understand that the coupling of phonon 
waves and electromagnetic waves is effective only for 
q — !> 0, since the speed of sound is negligible if compared 
with the speed of light. By combining a continuous ap- 
proximation for the description of the mechanical waves 
of the optical modes and Maxwell's equations for elec- 
tromagnetic waves one obtains the polariton dispersion 
curves. For a polar crystal the total dielectric function, 
which includes electronic and lattice polarization terms, 
then reads (here for simplicity in the scalar form): 



e{q-^ 0,w) = ec(g ^ 0,a;) + eiat(g ^ 0,a;) 



(2) 



where ec('Z -^ Ojw) is the long wavelength limit of the 
electronic contribution to the dielectric function (calcu- 
lated, e.g., in the RPA approximation), and eiat is defined 
as the contribution to the dielectric function due to lat- 
tice polarization. This latter quantity can be related to 



the electronic screening by a frequency-dependent gener- 
alization of the Lyddane-Sachs- Teller relation, that reads 



eiat(9 ^ 0,w) = ec(q -^ 0,w) 
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w|o - (w + ii]Y 



(3) 



In this last expression, the infinitesimal 77 -^ 0+. 

In the following discussion, and to keep the notation 
of the equations at a reasonably simple level, we will 
consider only cubic crystals with two atoms in the unit 
cell. However, all the equations can be easily generalized 
to arbitrary unit cells with more than two atoms, i.e. 
more than one infrared active modcj^ 

Equations ^ and ^ can be used directly within 
the GW framework without any further approxima- 
tion. However, the very large difference between typi- 
cal phonon (of the order of the meV) and electronic (of 
the order of the eV) frequencies makes the use of meth- 
ods like the contour deformatioii^i unpractical. We will 
therefore resort to the popular plasmon-polc model^i^ of 
the GW theory in order to derive a practical framework 
to include the effects of lattice polarization. 

For convenience, we now move to the reciprocal lattice 
G space and write e as a G, G' matrix. The electronic 
plasmon-pole model for e^^^Q, (w, </) reads: 



^cGG 



(q^^) 



JGG' 
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j'^ - [^GG'{q)-ivY 



(4) 



The two parameters J^qq/ (q) and lucg' (q) can be fixed 
for example by calculating £^(j(ji ('Z, w) at w = and 
Lu — iwp,— where wp is the plasmon frequency. In order 
to determine the total screening, in the case of q — > 0, 
one has to evaluate e^ G' (? ^ ^t^) by inverting the sum 
of the two matrices ec{q — > 0,aj) + eiat(q — > 0,0;). The 
way to proceed is to apply the matrix equality^^ 
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where 



(G + H)-^ ^ G-^ G-^HG 



g = Tr{HG-'} 



(5) 



(6) 



and G is an invertible matrix. In our case the matrix 
G is ecG,G'{q -^ 0,w) and H is eiat G,G'(g -> 0,uj). A 
useful simplification comes from the fact that H has only 
one non-zero element, the head of the matrix eiat oo{q -^ 
0,0;). In the following we will remove the explicit depen- 
dence on q of the dielectric matrix elements as it is clear 
that the lattice polarization contribution exists only for 
q — > 0. Using ^ in ([5]) it is easy to show that 



^GG 



G'(^) = e^GG'(^) 



Go(^)ecdG'(^) ^L 



SGO 



Urj,Q 



m('^) ^lo^i^+^vy 



(7) 



If we now replace e^ ^q, (w) by the standard plasmon- 
polc model (|H))^ we obtain our final expression for 



TABLE I. Calculated band gaps (in eV) for the selected crys- 
tals using different GW schemes compared with experimen- 
tal (exp) photoemission gaps. The experimental values of 
ojlo and uto (in mHartree) are also given. The experi- 
mental references are given in the first column. The labels 
GW'"-* and CW'^^'OLDA are abbreviations for scGVF'''' and 
GW^'^''@QPscGW , respectively. The last two rows represent 
the mean absolute error (MAE) and mean absolute percent- 
age error (MAPE). 

a;LO WTO GW@LDA GW'''* QPscGW scGW"'"-* 



The matrix elements of the wings have very similar 
expressions, thanks to the simplification of e^oo('^) ^P" 
pearing at the denominator in ([7|: 



exp 



LiF^ 


2.99 1.39 


13.24 


12.05 


15.81 


13.69 


14.20 


LiCl2^ 


1.74 0.87 


8.60 


7.97 


10.28 


9.05 


9.4 


NaCl2i 


1.21 0.75 


7.73 


7.14 


9.52 


8.37 


8.5 


Mg02£ 


3.29 1.82 


6.97 


6.38 


8.94 


7.71 


7.7 


Aip2a 


2.28 2.01 


2.32 


2.26 


2.79 


2.70 


2.49 


AlAs^ 


1.84 1.65 


1.88 


1.83 


2.34 


2.26 


2.23 


GaAs^ 


1.30 1.22 


1.16 


1.15 


1.52 


1.46 


1.52 


MAE 




0.60 


1.04 


0.74 


0.18 




MAPE 




11.5% 


16.3% 


9.4% 


3.0% 
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'^LO - ('^ + W)^ 



(10) 



with an analogous expression for the OG terms. 

The other matrix elements of ^qqi (w) are more compli- 
cated: beside the extra poles at w = ioiLOj there appear 
terms involving poles at uj = ±LJco, ^ = ±'^og, and 

LO = ±CJoo/eoo. 

We can then insert @ and ([7]) in the expression for 
the screened Coulomb interaction W = e~^v, where v 
is the bare Coulomb interaction. This allows to evalu- 
ate analytically the convolution integral of the self-energy 



Y.{q,uj)^— I doj'e"^^ G{q, uj + u;')W{q, uj') , (11) 



^GG'i^)- "^^^ resulting formulas are quite complex, but 
can be simplified for the head and the wings of the inverse 
dielectric matrix, remembering that phonon frequencies 
arc much smaller than the frequencies of electronic exci- 
tations. The head of the inverse dielectric matrix reads: 



eooH = l + - 

i. 

where we used 



^'00 



■'LO '^TO 



(woo - ivf ^°° '^LO - (w + ir?)2 ' 
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remembering that 



Gir,r',uj) 



E 



^kArWuAr') 



(12) 



w - ek,j + ir] sign(£fc J- - ^) 

where (j)k,j are Kohn-Sham or quasi-particle orbitals, ek,j 
are the corresponding energies, and fi is the Fermi level. 
The frequency integration (jll[) is performed using the 
residue theorem, in strict analogy with the procedure fol- 
lowed for the electronic screening alone, with the only 
difference that more than one pole is present in the gener- 
alized plasmon-polc model that we derived. Once again, 
the full result is rather cumbersome, but a relatively sim- 
ple expression is obtained for the lattice contribution to 
the matrix elements of S when only the head of the ma- 
0, uj) is retained: 



trix Er^r^, {q 
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')k',j\^liit\(Pk',j) 



AttujIo 



V 



IVrpQ 



E' 



0{fi-ek',j) 



^fc'jJ 



. ^^ hm 
(wLo - W 9-*-° 



P*k'.^AQ^O)pk'aJ{q^O) 



(13) 



with V the volume of the unit-ceh and pk',ij[q ^' 0) = 
limq^o/.!^fir-0*,_^(r)exp(-''3''')0fc/+gj(r). Note that 
this term only contributes to the occupied states due to 
the presence of the Heaviside function 0(/i — £k',j)- 

We implemented in ABINIT— the complete contribu- 
tion of the lattice polarization to E, i.e., the contribution 
of the head, wings and body of the dielectric matrix. 
However, only the poles at w = ±wlo can contribute 
significantly to the integral, as the coefficient in front of 
the other poles is always extremely small. For all sys- 
tems that we studied, the contribution of the head of the 
matrix alone is enough to describe the correction to the 
band structures within less than 1 meV, but this may not 
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be the case for non-cubic systems. 

Our calculations start by the determination of the 
standard Kohn-Sham ground-state, using the LDA and 
norm-conserving pseudopotentials. The LDA Kohn- 
Sham states are then used to calculate GW^QLDA band 
structures and as a starting point for QPscGVF runsi^ 
There are many different versions of GW, and conse- 
quently there are many different ways of using Eiat in the 
context of GW. We chose as starting point a converged 
QPscGW calculation, as this is in our opinion the most 
accurate level of theory available in ABINIT. However, 
one can not use Eiat in a self-consistent scheme, as this 
would induce contributions of the lattice polarization in 
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FIG. 1. (Color online.) Band structure of LiF along the 
highest symmetry lines of the Brillouin zone. From left to 
right: Kohn-Sham LDA bands, GW@LDA bands, "standard" 
QPscGVy bands, and GW'="@QPscGW bands, including lat- 
tice polarization in the screening. The direct band gap value 
is also given. The experimental gap is indicated with a hori- 
zontal red line. 



Eo, in disagreement with Eq. ^. We therefore decided 
to apply the corrections due to the phonons in a final 
perturbative step, in what we call the GW^^*^@QPacGW 
method. 

In the following, we present calculations for the band 
structures of a set of binary compounds: LiF, LiCl, NaCl, 
MgO, AlP, AlAs, GaAs. These are all polar sp materi- 
als, for which perturbative GW gives already good band 
gaps, with varied bond ionicities, and band gaps ranging 
from 1.52 eV to 14.2 eV. All calculations were performed 
with ABiNlTi^ The energy cutoff was set to 50 Ha for 
LiF, LiCl, NaCl, and MgO, and to 30 Ha for the remain- 
ing compounds, and the fc-point grids were a 6 x 6 x 6 
Monkhorst-Pack. This allowed for a convergence in the 
total energy to better than 1 mHartree. The cutoff for the 
dielectric constant was set to 5Hartree, and 142 empty 
states were used to obtain a convergence in the energy 
gaps to better than 0.1 eV. Finally, and in order to limit 
the effect of other possible sources of imprecision, we de- 
cided to use for these calculations experimental lattice 
constants and experimental LO and TO phonon frequen- 
cies at r (see TableU]). Of course the method can be easily 
made fully ab initio by calculating such quantities within 
DFT,^ 

In Table |T] we report the calculated band gaps for 
the selected binary compounds using the different GW 
schemes. As it is well known, GM^@LDA underesti- 
mates the band gaps of sp compounds by slightly more 
than 10% while QPscGW overestimates them by about 
10%!^ Turning on the correction due to lattice polariza- 
tion cuts down the error of these techniques by a factor 
of three, bringing the mean absolute percentage error to 



3%. These results clearly prove that the overcstimation 
of the band gap by QPscGM^ is not mainly due to a 
deficient treatment of electronic correlation, but to a de- 
ficient description of the screening of the medium. We 
remark that the inclusion of the lattice contribution to 
the screening in a G1^@LDA calculation deteriorates the 
agreement with experiments, leading to a mean absolute 
percentage error of 16%. This shows that the good per- 
formance of perturbative GW calculations for sp polar 
materials is due to a partial cancellation of errors: the 
underestimation of the band gap opening is compensated 
by the neglect of the band gap shrinkage due to the lat- 
tice polarization. 

In Fig. [1] we display the band structures of LiF ob- 
tained with the different approximations under study. 
We can observe that the inclusion of the lattice polar- 
ization in the screening produces a rigid shift downwards 
of the conduction bands, with negligible effects on the 
band widths and band dispersions. While the same qual- 
itative behavior is observed for all the materials under 
study here, we cannot exclude that other effects may be 
present in more complex materials. 

Finally, we would like to make two remarks: (i) The 
coupling between electrons and phonons in our model 
is indirect and comes through Maxwell equations. It is 
therefore completely unrelated to electron-phonon cou- 
pling, (ii) To take into account all phonon contributions 
to the electronic band structure, both the lattice contri- 
bution to the screening and the Fan and Debye- Waller 
(and possibly higher order) terms of the AUen-Heine- 
Cardona theorj*i^ should be included. It is true that for 
the materials studied here one docs not expect a large 
contribution from the electron-phonon coupling, but it is 
not inconceivable to find a strong ionic compound with 
strong electron-phonon coupling. Note that both contri- 
butions will tend to decrease the purely electronic gap. 

In conclusion, we developed a fully ab initio GW 
framework that includes the effects of the screening ow- 
ing to the polarization of the lattice. Within this frame- 
work we show that the overcstimation of the band gaps 
by restricted self-consistent GW techniques is due to the 
neglect of this contribution, and we manage to bring the 
error in the calculated band gaps to a mere 3%. The 
lattice contribution decreases the band gap by a factor 
that increases with the size of the longitudinal-transverse 
optical splitting and it can represent more than 15% of 
the band gap in highly polar compounds. These results 
call for a reexamination of many theoretical calculations 
of the quasi-particle spectrum for polar materials, includ- 
ing for many oxides, nitrides, etc. that are important in 
several fields of technology. 
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